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. Abstract. This review discusses (i) dynamical methods for determining 

I the masses of Galactic and extragalactic star clusters, (ii) dynamical pro- 

. cesses and their time-scales for the evolution of clusters, including evap- 

1-^ I oration, mass segregation, core collapse, tidal shocks, dynamical friction 

■ and merging. These processes lead to significant evolution of globular 

7—i [ cluster systems after their formation. 
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. !• Introduction 

(N ■ 

I The Milky Way and probably all large galaxies contain old globular cluster 

^ ' populations (see the review by Harris 1991). These old star clusters have an 

I approximately log-normal luminosity function, and the mean cluster luminosity 

is somewhat brighter than My = — 7 with little dependence on the host galaxy 
luminosity. In the Milky Way their typical mass-to-light ratios are M/Ly ~ 2, 
and typical total masses are ~ 2 x 10^ M0 (Pryor & Meylan 1993). It is widely 
2 ! assumed that the globular clusters we see today must be the part of an initially 

' larger population that survived the internal and external dynamical processes 

. leading to cluster destruction (e.g., Ostriker 1988). 

' One of the exciting results from HST has been the discovery of young star 

. clusters in starburst and interacting galaxies. Whitmore & Schweizer (1995) 

/\ ' found many hundreds of young clusters in the Antennae galaxies. Young cluster 

^ . systems have now been discovered in other interacting and merging galaxies, 

" " ' in barred and starburst galaxies, and even dwarf starburst galaxies (e.g., ESO 

338-IG04, Oestlin, Bergvall & Roennback 1998). The luminosity functions of the 
young clusters are not log-normal, but seem to be better described by power- 
laws, about cx L~^. Carlson et al. (1999) use population synthesis models to 
determine the ages of the blue clusters in the young merger remnant NGC 3597. 
Based on these models they argue that the difference in the observed luminosity 
function when compared to the Galactic globular clusters cannot simply be 
an age effect, even if the young clusters formed with an intrinsic age spread. 
Are these young cluster systems then a good model for what the Milky Way's 
globular cluster population could have looked like at birth? 

This review gives a brief discussion of dynamical methods to determine 
masses of distant and nearby star clusters (Section 2). It then goes on to de- 
scribe a number of dynamical processes and their time-scales which will lead 
to evolution and potentially destruction of star clusters over long time-scales. 
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Finally, the results of some evolution calculations for globular cluster systems 
are briefly summarized (Section 3). 

2. Dynamical Mass Determination for Star Clusters 

In this Section we discuss methods for estimating star cluster masses from struc- 
tural and kinematic measurements. Mass estimates based on stellar population 
properties are discussed in U. Fritze von Alvensleben's article in these proceed- 
ings. 

2.1. Virial Masses 

A simple global mass estimate for a star cluster can be obtained from the virial 
theorem. This says that, in equilibrium, the radius of a stellar system is propor- 
tional to GM/V"^, where M is the total mass and V the rms three-dimensional 
velocity of the stars. The constant of proportionality generally depends on the 
stellar density profile, but Spitzer (1969) showed that if the relation is expressed 
in terms of the half-mass radius r/i, this dependence is weak and the constant is 
approximately 0.4 for realistic cluster profiles. If we furthermore assume that the 
cluster is spherical, V"^ = Sfjp where a\\ is the one-dimensional rms velocity dis- 
persion along the line-of-sight, and write cJio = <T||/10kms~^ and = r/i/5pc, 
then 

Mv = ^.^a\rnlG = 8.7 x lOVf^rg M©. (1) 

When using this formula to estimate star cluster masses from observed velocity 
dispersions and radii, a few points should be noted: 

(i) Because the virial mass (||) is a global estimate, it is independent of ve- 
locity anisotropy. For example, shifting some stars to radial orbits while keeping 
the (spherical) potential fixed, will result in a larger central velocity dispersion 
but also lead to reduced velocities in the cluster halo. To maintain virial equi- 
librium these changes must add in just such a way that the global (T\\ remains 
the same. 

(ii) The dynamical evolution of star clusters leads to mass segregation and 
the formation of a halo of low-mass stars on preferentially radial orbits (see 
§ 3.2.| below). For evolved clusters the measured half-light radius will therefore 
in general underestimate the half-mass radius, the observed velocity dispersion 
will underestimate the rms velocity dispersion, and eq. (0) will underestimate 
the mass. 

(iii) Sometimes only the velocity dispersion for stars in the core is known, 
or the velocity dispersion from integrated light within some aperture. In these 
cases, a dynamical model is needed to convert this to the rms (T||. This intro- 
duces some uncertainty in the mass estimate because the derived (Th depends on 
anisotropy. 

With high-resolution spectra and HST photometry virial masses can be 
determined for some young 'superclusters' seen in starburst galaxies. Masses for 
three clusters in M82 are compared by Smith & Gallagher in these proceedings, 
spanning a range from 3 x 10^ Mq — 2 x 10^ M©. 
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2.2. Core Masses 

Rood et al. (1972) gave a formula that is often used to determine core masses. 
This is based on the dynamics of King models (King 1966; Binney & Tremaine 
1987) and assumes that the velocity distribution in the core is isotropic: 



^\ = (2) 

Here ctq is the central velocity dispersion, Iq the central surface brightness and 
the core radius. The product of the two last quantities is rather insensitive 
to errors caused by seeing. Eq. (P) is very accurate as long as the assumption of 
isotropy is met (Richstone & Tremaine 1986); but it can overestimate the mass 
by a factor ~ 2 if the system is actually radially anisotropic (Merritt 1988). 

Core mass-to-light ratios can only be determined for well-resolved Galactic 
star clusters for which the core parameters Tc, Iq and ctq can be estimated. Even 
with the resolving power of HST the cores of distant young clusters cannot be 
resolved. 



2.3. Masses from Model Fitting 

An alternative method of estimating cluster masses is fitting the photometric 
and kinematic data with dynamical models. A simple such scheme was used 
by Djorgovski et al. (1997) in their study of the M31 globular clusters mass-to- 
light ratios. They used structural and photometric parameters for these clusters 
obtained with HST and kinematic measurements in a rectangular aperture ob- 
tained with Keck and HIRES. They then estimated an aperture correction from 
King models to transform their measurements to central velocity dispersions, 
and used a formula analogous to eq. (^ to estimate masses with the constant 
again determined from models. 

For some Galactic globular clusters large velocity samples are available and 
in such cases much more detailed model fitting is possible (Pryor et al. 1989). 
The masses of the Galactic globular clusters referred to in §^ (Pryor & Meylan 
1993) have been determined by these techniques. A recent such study is Cote 
et al. (1995) who investigated the dynamics of the globular cluster NGC 3201, 
using a CCD surface brightness profile and a sample of 857 measured stellar 
radial velocities to trace the velocity dispersion profile to large radii. 

In such work the data are fitted by single- or multi-mass King-Michie mod- 
els. In the multi-mass models, a power-law mass function for the cluster stars 
is typically assumed, and for each mass bin mj, a distribution function of King- 
Michie type (Michie 1963, Binney & Tremaine 1987) is used: 

/,(i?,J)oce-^-^'(e-^»^-l). (3) 

Here E, J are the specific energy and angular momentum of a star in the (spher- 
ical) star cluster potential, and /3 can be thought of as specifying an anisotropy 
radius. In the core of the cluster, stars of different masses are assumed to be in 
equipartition (§ ^.3j ), so that Ai (x rrii. In the fitting procedure the free param- 
eters are the radius, velocity and luminosity scale, the cluster's concentration 
parameter, the anisotropy radius, and the index of the mass function. These 
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parameters are determined from fitting to the measured surface brightness and 
velocity dispersion profile. This leads to a determination of the mass-to-light 
ratio profile and anisotropy profile, rather than just a single M/L as for the 
previously described techniques. However, the fit is non- unique in the sense 
that adding even fairly large numbers of faint low-mass stars in the halo (ex- 



pected there from mass-segregation and evaporation, see §3.2. below) have little 
effect on the observed profiles. In their study. Cote et al. (1995) find a steady 
rise in M/L with distance from the cluster center, as expected from dynamical 
evolution theory, and a global M/Lb — M/Ly — 2.0 it 0.2. 

2.4. Masses from Proper Motions 

For nearby Galactic globular clusters, it is possible to measure stellar proper 
motions in addition to radial velocities. Proper motion measurements give in- 
formation about the velocity dispersions in two directions on the sky (radial 
along projected R, and tangential), and for a spherically symmetric cluster they 
are therefore in principle sufficient to determine the velocity ellipsoid as a func- 
tion of radius, and thus the mass profile free of assumptions about anisotropy. 
The projected proper motion dispersions aji{R) and axiR) are related to the 
intrinsic velocity dispersions ar{r) and (Tt{r) by Abel integral equations and 
can thus be inverted (Leonard &; Merritt 1989). Moreover, from the inferred 
(Jr{r) and iTf(r) one can predict the line-of-sight velocity dispersions C||(i?) and 
compare with independent radial velocity data. This provides a check on the 
modelling and also can be used to determine the cluster distance. In terms of 
global velocity dispersions, (cr^) = ((o"jj) + {aj))/2 for a spherical cluster and 
the correct distance. 

In an early study along these lines Leonard et al. (1992) investigated radial 
velocity and proper motion data for the globular cluster M13. They concluded 
that the mean anisotropy of this cluster (3 = 3((T^ — o"|n)/(3<T|j — (t|.) ~ 0.3 and 
that the effect of the anisotropy on the mass determination is ~ 20%. Much 
more detailed modelling will be possible with the large proper motion surveys 
currently in progress. 

2.5. Non-Parametric Cluster Mass Distributions 

With large samples of stellar velocities at hand, radial velocities or proper mo- 
tions, it is possible to infer the mass distribution of the cluster without making 
specific assumptions like King-Michie stellar distribution functions. This re- 
quires solving the Jeans and projection equations for the intrinsic density and 
velocity dispersions under some smoothness constraint, given the data. I do not 
give the equations here, but refer to the papers mentioned below. 

When the data consist of several hundred radial velocities, some assumption 
about the anisotropy is still needed. Gebhardt & Fisher (1995) describe such a 
non-parametric analysis of radial velocity data for four Galactic globular clusters, 
assuming isotropy of the stellar orbits. With a few hundred stellar velocities in 
each case the results are still noisy, but indicate radially increasing M/L-profiles 
as expected. Merritt, Meylan & Mayor (1997) describe a similar analysis of the 
cluster u) Centauri, assuming that it is oblate and seen edge-on, and that it is 
described by a meridionally isotropic two-integral model. They find that the 
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mass distribution cannot be strongly constrained by their data, but appears to 
be slightly more extended than the luminosity distribution. 

As discussed above, proper motion data result in two independent velocity 
dispersions in the plane of the sky, and thus, within a spherical model, they con- 
tain sufficient information to determine the anisotropy of the stellar orbits. With 
sufficiently large data sets it will therefore be possible to model the anisotropy 
profile and mass distribution of a spherical cluster non-parametrically. 



3. Dynamical Evolution Processes and their Time-Scales 
3.1. Relaxation 

On dynamical time-scales large star clusters {N >> 100) evolve collisionlessly. 
That is, for some time after birth they are described by a quasi-equilibrium 
phase-space distribution function (df) which is a function of the integrals of 
motion (or of the stellar orbits) in the mean field potential. On longer times- 
scales, however, the graininess of the distribution of stars becomes important, 
and the dynamical evolution is no longer collisionless. Over a relaxation time 
two-particle interactions then deflect the cluster stars from the orbits they would 
otherwise have followed in the mean gravitational potential. 

In the approximation of a homogeneous distribution of equal-mass stars, 
with density pQ and isotropic Maxwellian velocity distribution with dispersion 
ctQ) the two-body relaxation time is (Spitzer & Hart 1971) 

= n2 1 A = ^ ('wn.'m-l {lnA)i,^ yr. (4) 

Lr nijf po In A 

Here n = 10^n4pc~^ is the number density of stars, ctq = lOcriokms"^ the 
velocity dispersion, m,, = m^,^Q Mq is the mean mass per star, A is of order the 
number of stars, and (In A)io = (In A)/10. Note that trfi is inversely proportional 
to the stellar phase space density. Central relaxation times in globular cluster 
cores evaluated with eq. (Q) are ~ 10^ — 10^ yr. 

The relaxation time often varies by large factors between the central and 
outer parts of a stellar system. It is then useful to define a half-mass relaxation 
time. For a virialized star cluster this is obtained from eq. (§) by replacing p^ 
with the mean density inside the system's half-mass radius = Srs pc and ctq 
by one third of the rms V'^ , and then using the virial theorem to express V 
through r/j and the total cluster mass M = lO^Ms Mq. The result is (Spitzer & 
Hart 1971) 

0.14iV ( rl \ ^ N s 1/2 3/2 1 . ^ 1 

t^^ = _ -A- = — _ td = 7.2 X 10^ M;^^ r'^ In A yr, 

ln0.4iV \GMj 26 1og0.4iV 5 5 *,© v /lo J > 

(5) 

where is the total number of stars in the system and 

U = Th/V = 1.58 (rl/GM)^/^ = 8.3 x lOVg^/'Mg-'/' yr (6) 

is the dynamical time. For comparison with the local formula eq. (^, the fiducial 
values used in eq. (|5|) correspond to a one-dimensional virial velocity dispersion 
a ~ 3.4kms~^ and a mean density n ~ 96pc~^. 
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On short time-scales cluster evolution is still collisionless, so long as td <C trh 
(requiring N » 100); for example, during the violent relaxation at formation. 
The resulting quasi-equilibrium df subsequently evolves slowly in response to 
collisions, which will tend to drive the system towards an isothermal energy 
distribution. One aspect of such slow evolution would be a decrease of ellipticity 
with dynamical age (Fall &: Prenk 1985). This could be the reason for the 
significantly rounder globular clusters in M31 and the Milky Way as compared 
with the LMC and SMC clusters (Han & Ryden 1994). 

3.2. Evaporation and Core Collapse 

Collisions between single stars modify the stellar df in two ways. The rarer 
process is ejection, in which a single close encounter leads one of the stars to 
acquire a velocity greater than the local escape velocity Vg and to escape from 
the cluster. The time-scale for this is tej = -N/{dN/dt) ~ 1.1 x lO^lnOANUh 
~ W^trh (Henon 1969). The more important process of evaporation is caused 
by the cumulative effect of many weak encounters, which gradually increase a 
star's energy until v > Vg. It is easy to show that the rms escape velocity of 
the cluster is just twice its rms virial velocity. Thus, on average, a particle with 
V > IV = \/l2c7|| will escape. For a Maxwellian velocity distribution, a fraction 
e ~ 0.74% of stars have v > 2V; these stars will escape in one dynamical time, 
after which the high- velocity tail is repopulated only in ~ trh- Thus one expects 
the evaporation time scale of the cluster to be ~ {e/trh)~^', detailed calculations 
(Spitzer & Thuan 1972) show that 

tev = -N{dN/dt)-^ ~ 300 trh- (7) 

Because the evaporation is dominated by weak encounters, escaping stars leave 
the cluster with only very small positive energy; thus the total energy of the 
remaining cluster is nearly constant, but must be shared among a shrinking 
number of stars. In virial equilibrium N'^/rh — const, and thus p oc N/rf^ cx 

and trh oc Nrl^^/M^/^ oc N"^/^. So as the cluster becomes denser, evaporation 
accelerates and the system contracts to negligible mass and radius in finite time. 

This evaporation model, however, neglects the fact that the evolution of the 
stellar cluster is not homologous and that the rate of evolution is much faster in 
the dense core than in the system's outer parts. Stars gaining energy towards 
evaporation build up an extended halo where the time scale for further energy 
gain increases strongly, so that these stars may not in fact escape during the 
age of the cluster. On the other hand, the dense core loses stars to the halo on 
the much faster central relaxation time, and may collapse to very high densities 
before Mtot and rh can change much. 

This phenomenon of core collapse may be understood as a consequence of 
the fact that self-gravitating star clusters have negative specific heat (Lynden- 
Bell &: Wood 1968): In virial equilibrium the total energy E = —T, where 
T is the total kinetic energy, which is proportional to the virial temperature 
T/M = V^. As energy is withdrawn from the cluster, its kinetic energy increases 
and so does the virial temperature. Since rh — GM'^/(2\E\), the cluster thereby 
contracts. Vice- versa, an energy production mechanism (e.g., from binary stars) 
causes the cluster to cool and expand. Now the dense core of the cluster may be 
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approximately regarded as a virialized system in thermal contact with the rest 
of the cluster. It is normally hotter than its surroundings and therefore loses 
energy to them through stellar encounters. As a result it shrinks and becomes 
yet hotter, loses still more energy to the surrounding stars, and contracts to 
formally zero radius in finite time. 

For a single-mass star cluster the late stages of core collapse are self-similar 
(Lynden-Beh & Eggleton 1980, Cohn 1980). As the core radius 



Tc = SaJ^AnGpc (8) 

shrinks, the central density pc and velocity dispersion cjc increase and the core 
mass Mc decreases according to 



fje oc oc r-o-ii, M,(x pcT^^r'^/'^ (9) 

until the core radius and mass formally shrink to zero at time tec- Moreover, the 
density profile of the cluster outside the collapsing core has the same exponent: 
p oc r~^-^^ for rc{t) <C r <C rc(0). 

The time-scale for core collapse is proportional to the central relaxation 
time; for a single mass cluster it is tec — 330 t^c once the collapse is in the 
self-similar phase (Cohn 1980, Heggie &: Stevenson 1988). The total time until 
core collapse in Cohn's (1980) model is ~ 16 half-mass relaxation times or ~ 60 
initial trc- As Goodman (1993) has emphasized, the former number depends 
on the mass distribution of the cluster, and tcc/trh will be less than 16 for 
clusters more centrally concentrated than Plummer models. By noting that 
r~^drc/dt oc t~^ oc Pc/cTc r~^'^^, one can solve for the asymptotic time- 
dependence of the collapse: 

oc (t,e - ^)°■^^ Pc ^ {tec - tr^-^\ cTcoc (^cc-^)"°■°^ Mc ^ {tec - tf^\ 

(10) 

In summary, the collapse of a single mass cluster occurs in two stages (Cohn 
1980). The longer part of the evolution is an evaporative phase, during which 
stellar collisions simultaneously populate a halo and make the core shrink and 
become denser. Only towards the end does the evolution accelerate and enter 
the gravothermal instability phase of self-similar collapse. 

3.3. Equipartition, Mass Segregation, and Multi-Mass Core Collapse 

When the cluster contains different stellar masses, energy can flow not only from 
the core to the halo, but also between stars of different masses. Stellar collisions 
drive the system towards equipartition of energy mi{vf) = const. As eq. (^ 
shows, relaxation proceeds faster for more massive stars. The equipartition 
time-scale measures the rate at which a group of heavy stars with masses m2 
loses energy to lighter stars of mass mi (Spitzer 1969): 

teq - Q/^ ~ ^'^ tMmi) (11) 

8(D7r)^/^G^mim2ni m A m2 

where we have assumed equal temperatures (vf) = (ul) and used eq. (^. Ini- 
tially, (vf) is independent of stellar mass; thus the massive stars lose kinetic 
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energy and sink to the center, while hghter stars gain kinetic energy in colhsions 
and move outwards, a process caUed mass segregation. Moreover, eq. ([l^) shows 
that mass segregation of the massive stars occurs before relaxation of the cluster 
as a whole becomes significant. 

However, equipartition may never be reached. A simple case considered by 
Spitzer (1969) is one with two mass groups such that the heavy stars are much 
more massive than the light stars, m2 ^ mi, but the total mass in the cluster 
core is dominated by the light stars: M2 <C Pir^i- In this case equipartition 
causes the formation of a small subsystem of heavy stars (M2, 1112) in the core of 
the distribution of lighter stars. Applying the virial theorem to the subsystem 
of heavy stars gives 

a4GM,^4^ (12) 

where the first term describes the self-energy of the subsystem M2 and the second 
term its interaction with the system of light stars (fc is a constant of order unity). 
Spitzer (1969) noticed that the right-hand-side of eq. (^) has a minimum when 
regarded as a function of rh2- An equilibrium can therefore exist only if (w^) is 
greater than this minimum, that is, assuming equipartition, if 

J% < 4.0.- (^)"\ (13) 

In other words, if its mass is too large, the subsystem of heavy stars remains a 
dynamically independent stellar system with mean square velocity greater than 
the equipartition value. It continues to lose energy to the lighter stars, becoming 
denser and hotter, and evolving away from equipartition all the time {mass 
stratification instability). Fokker-Planck calculations (Inagaki & Wiyanto 1984, 
Cohn 1985) show that in the end the subsystem of heavy stars core collapses 
independently from the cluster of light particles, just like a single mass system. 

The evolution to core collapse with a spectrum of stellar masses has been 
considered by Inagaki & Saslaw (1985) and Chernoff & Weinberg (1990). The 
detailed evolution occurs in several phases: First, collisions trying to establish 
equipartition of energy lead to mass segregation and the formation of a heavy 
mass core. Then this core undergoes the gravothermal instability, i.e., contracts 
while remaining hotter than the mean temperature of the system and conducting 
energy outwards. This collapse accelerates towards core collapse, and finally 
goes over into a single-component collapse which reaches formally infinite central 
density. The time scale for this multi-mass core collapse evolution is faster than 
that for core collapse in any single component cluster, typically a factor of a few 
faster than for a cluster composed of the heaviest mass alone. 

Deep in collapse, the density slopes of all mass groups m^ are characterized 
by separate power laws in the region where the heaviest component dominates 
the potential, such that approximately (Cohn 1985, Chernoff & Weinberg 1990) 

d\npk/d\nr ~ -1.89(mfe/m„) - 0.35, (14) 

where m^ is the mass of the heaviest species in the cluster. The overall mass 
profile is then not self-similar. 
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A multi-mass core collapse, however, may be strongly influenced by the 
stellar evolution of the more massive stars. This has two main effects: First, 
the mass loss from massive stars through winds may lead to an overall mass loss 
from the cluster, and thus cause an adiabatic expansion. Secondly, the finite 
stellar life-time tMS limits the time tec during which they can core collapse, such 
that tcc{m^) < tMs{iT^*)- Both effects greatly increase the overall core collapse 
times; compared to a system of point masses within the range (0.4 — 15)m0, 
Weinberg & Chcrnoff (1989) find an increase by about a factor of 30 — 60 in 
their globular cluster models, including the expansion effects. 

A reasonable approximation for the stellar lifetime of all but the most mas- 
sive stars is tjv/5 — 9 • 10^(m*/mQ)~^'^ yr. Thus if the most massive stars leave 
black hole remnants of SM©, these together with tight binaries will dominate 
the evolution after 5 ■ 10^ yr, while if the most massive remnants are 1.4 
neutron stars, they and the binaries will dominate after 4 • 10^ yr. The latter 
time scale approaches the time expected for core collapse in typical Milky Way 
globular s. 

3.4. Reversing core collapse 

A number of energy source mechanisms can stop core collapse (e.g., Goodman 
1993): (i) Processes that generate kinetic energy in the core directly, such as 
binary stars transferring energy to the field stars in collisions, (ii) Mass loss 
processes that heat the core indirectly above its virial temperature, including: 
normal stellar evolution, accelerated stellar evolution by the formation of massive 
stars in mergers, and ejection of stars through binaries. In all cases the net result 
is adiabatic expansion and cooling of the core. 

Only hard binaries contribute to field star heating. Binaries are hard if 
their binding energy = —Gmim,2/ a (with a their semi-major axis) exceeds 
the mean kinetic energy, E^, > Sm^a^; those with Ef, < Sm^a^ are called soft 
binaries. Reggie's law (Reggie 1975, Rut 1983) states that, on average, hard 
binaries get harder by collisions with field stars, and soft binaries get softer. 
Essentially, the orbital velocity of a hard binary is on average greater than the 
velocity of an incoming field star, and the tendency towards energy equipartition 
therefore results in a net transfer of energy to the field star. The opposite is 
true for soft binaries, which gain net energy and eventually dissolve. The binary 
behaves like a mini-system with negative specific heat: as energy is withdrawn 
from it, the orbit shrinks, the orbital velocity increases, and the binary hardens. 
When the binary becomes sufficiently hard, the typical recoil from a collision 
with a single star becomes large, and the binary will eventually be kicked out 
of the cluster. Just like in the Sun, the binaries providing the nuclear energy 
source will eventually be 'burned'. 

Binaries can be formed by a close gravitational interaction of three stars 
('three-body binaries'), by dissipational tidal capture, or at the time of star 
formation ('primordial binaries'). To be effective in reversing core collapse, 
binaries must have orbital semi-major axes 

a<Gm*/3a^ = 3 aj^^ m*^© AU. (15) 

The formation of a hard three-body binary requires a close encounter be- 
tween two stars {5v ~ v) with a third star in the immediate vicinity, such that 
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one of the three stars acquires additional energy, leaving the other two as a 
bound pair. Thus the time-scale is (Goodman 1984, Binney & Tremaine 1987) 

ts (np^v)-^ (np^)-^ ~ ^— ~ InN^Uo, (16) 

where p ~ Gm^/v^, v = 0{a) because low relative velocities dominate, and we 
have used eq. (^) to express Goodman's result in terms of the central relaxation 
time and the total number of stars in the core, Nc- This implies that about 
1/Nc In A^c three-body binaries form per central relaxation time. In other words, 
three-body binaries become important if the final core collapse is driven by fewer 
than 100 of the largest mass stars. 

3.5. Tidal field and tidal shocks 

A steady tidal field lowers the energy threshold beyond which stars are no longer 
bound to the cluster. It thus increases the mass loss rates from evaporation, 
both because the fraction of stars in the velocity distribution that escape in a 
dynamical time increases, and because the decreasing number of stars in the 
cluster leads to shorter relaxation times. The Quintuplet and Arches young 
clusters (Figer et al. 1999) in the inner Galactic bulge are two clusters for which 
these tidal effects are very important (Kim et al. 1999). 

In reality, the tidal field is not stationary in the frame of the cluster stars. 
This complicates the escape process, but more importantly it leads to a new dy- 
namical process in cluster evolution, referred to as gravitational shocking (Os- 
triker, Spitzer & Chevalier 1972). The tidal field acting on the cluster may 
suddenly increase in strength when the cluster passes through the disk of its 
host galaxy, or when it comes close to the high-density inner bulge near peri- 
galacticon of its orbit. In both cases, the perburbations to the stellar orbits 
caused by the tidal shock lead to an effective energy input in the cluster which 
makes the cluster less bound and accelerates mass loss from internal processes. 

A detailed recent discussion of this process is given by Kundic & Ostriker 
(1995) and Gnedin, Lee & Ostriker (1999). For stars in the outer parts of 
the cluster, the tidal perturbation can be approximated as impulsive because 
of the short time-scale of passage through the galactic disk. In the cluster's 
central parts, on the other hand, the stellar orbital time-scales are short and 
adiabatic invariance reduces the effects of the perturbation. Traditionally these 
effects of the tidal shock were described by a first-order term {{/S.E)ts)^ which 
denotes the net energy gain averaged over stellar orbits at a given position in the 
cluster. Kundic &: Ostriker (1995) noticed that the second-order term {{AE)^^) 
is typically even more important and competes with two-body relaxation near 
the half-mass radius in driving evolution of the cluster's internal structure. This 
can speed up core collapse by a factor of three (Gnedin, Lee &: Ostriker 1999). 
Cluster destruction is accelerated; recent modelling of the evolution of the Milky 
Way's globular cluster system shows that the typical time to destruction becomes 
comparable to the typical age of the Galactic globulars (Gnedin & Ostriker 
1997). 
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3.6. Dynamical friction and merging 

As already noted by Tremaine, Ostriker & Spitzer (1975), massive star clusters 
experience dynamical friction against field stars as they move along their orbits 
through the host galaxy Because of the frictional drag the cluster loses orbital 
energy and spirals into the galaxy center, where the tidal field becomes ever 
stronger and will eventually dissolve the cluster. 

The time-scale for dynamical friction for a cluster on a circular orbit at 
initial radius rj = 2rj^2 kpc in a singular isothermal sphere with circular velocity 
Vc = 250^250 knis~^ is (Chandrasekhar 1943, Binney & Tremaine 1987) 

1 1 7r^7; 

where is again the cluster mass in units of 10^ Mq. The friction time-scale 
thus scales with the square of the cluster's initial radius in the potential, and is 
inversely proportional to its mass. It is the inner, most massive clusters which 
are affected first. 

If we continue to model the inner parts of the host galaxy as an isothermal 
sphere with Mg{t) = vIr/G and use the virial theorem [eq. (||)] for the cluster, 
we can determine the radius at which the incoming cluster will dissolve as 

1 



^dis = Th \^ — j = ^jY^^ = 46 rs cJio W250 pc. (18) 

Young clusters formed in the high-density regions of starburst galaxies would 
thus contribute to the build-up of the nuclear bulge after being dragged inwards 
by dynamical friction and tidally shredded by the tidal field. 

In some circumstances it may be possible that several young clusters are 
born close enough to eachother to tidally interact and even merge. To quantify 
this we use an approximate merging criterion fitted by Aarseth & Fall (1980) 
to the results of N-body merger simulations. For the escape velocity of the 
clusters at pericenter p of their relative orbit we take an approximate expression 
assuming two overlapping Plummer spheres, v1{p) = 27.6crj^/(l -|- p^/1.2r^)-'^/^ 
(see also the discussion in Gerhard & Fall 1983). Then the criterion for merging 
becomes 

where Vp is the relative velocity at pericenter. Here we have used the virial 
equation (|^), and cJ|| is again the one-dimensional rms velocity dispersion of the 
cluster. For head-on collisions, this formula predicts merging for Vp ^ 6cJ|| = 
60cJiokms~^ (slightly more than ^/2 times the rms escape velocity from each 
cluster), or Aw = \/36 — 27.6 (T|| ~ 3cr|| ~ 30o"iokms~^ for their relative velocity 
at large separations. It also shows that merging requires the two clusters to 
come within several half- mass radii of eachother for merging to occur, at cor- 
respondingly smaller approach velocities. The most likely situation for this to 
happen would be when two clusters are born from the same giant molecular 
cloud complex. 
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3.7. Evolution of globular cluster systems 

The evolution of globular cluster systems has recently been modelled in a num- 
ber of studies, among others by Gnedin & Ostriker (1997), Murali &; Weinberg 
(1997), Baumgardt (1998) and Vesperini (1998). These models combine as- 
sumptions about the initial cluster mass function and cluster locations with 
evolutionary models for individual clusters. In the models the various processes 
described above are considered, and treated in some studies by parametrized 
mass loss rates or analytic approximations to the results of N-body simulations, 
in others as diffusion terms in Fokker-Planck models. Some of the main results 
of these studies are: 

(i) Globular cluster systems in galaxies evolve significantly. In the Milky 
Way the typical cluster destruction time is of order the age of the system, and 
about half of the present globulars will be destroyed in the next Hubble time. 

(ii) Clusters in the inner regions of their host galaxy are disrupted most 
rapidly. Similarly, clusters on eccentric orbits are preferentially destroyed over 
clusters on tangential orbits. 

(iii) Low-mass and high-concentration clusters are disrupted by evaporation, 
loosely bound clusters and those on central or eccentric orbits by tides, and 
massive inner clusters by dynamical friction and tides. 

(iv) Low-mass clusters are destroyed most efficiently and initial power- 
law mass distributions tend to become transformed towards approximately log- 
normal mass distributions. 
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